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Abstract. The Jacobi-Stirling numbers were discovered as a result of a problem involving the spec- 
tral theory of powers of the classical second-order Jacobi differential expression. Specifically, these 
numbers are the coefficients of integral composite powers of the Jacobi expression in Lagrangian 
symmetric form. Quite remarkably, they share many properties with the classical Stirling numbers 
of the second kind which, as shown in [14] . are the coefficients of integral powers of the Laguerre 
differential expression. In this paper, we establish several properties of the Jacobi-Stirling numbers 
and its companions including combinatorial interpretations thereby extending and supplementing 
known contributions to the literature [2], [3], [7J, [II], and |16| . 



1. Introduction 
The Jacobi-Stirling numbers, defined for n,j € N by 

n\ ._v-, , Nr+i (a + P + 2r + l)r(q + (3 + r + 1) [rjr + a + (3 + l)] 7 



(L1) UJa/ S ( l ' r l( j -r)T(a + f3 + j+r + 2) 

were discovered in 2007 in the course of the left-definite operator-theoretic study [10] of the classical 
second-order Jacobi differential expression 

(1.2) t a>fi [y](x) = Xy(x) (a; €(-1,1)), 
where 

(1.3) iaAvM ■= ~T~: ( ("(I " + xf +1 )y'(x))' + k(l - x) a (l + xfy(x) 

= -(1 - x 2 )y"{x) + (a-/3 + (a + /3 + 2)x)y'(x) + ky{x) {x G (-1, 1)). 

Here, we assume that a,/3 > —1, A; is a fixed, non-negative constant, and w a ^{x) is the Jacobi 
weight function defined by 

(1.4) w a ,p(x):=(l-x) a (l + xf (i£(-l,l)). 

The Jacobi-Stirling numbers are the main subject of this paper; as we will see, these numbers 
are similar in many ways to the classical Stirling numbers of the second kind. Moreover, this 
manuscript may be viewed as a continuation of the combinatorial results obtained in [2], [3], and 
[7]; each of these papers deals exclusively with various properties of the Legendre-Stirling numbers 
{PSn^}, a special case of the Jacobi-Stirling numbers. Indeed, by definition, PSn^= {j} 00 - The 
Jacobi-Stirling numbers have generated a significant amount of interest from other researchers in 
combinatorics. In this respect, we note that the Legendre-Stirling numbers appear in some recent 
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work [4J related to the Boolean number of a Ferrers graph; these authors also show that there is an 
interesting connection between Legendre-Stirling numbers and the Genocchi numbers of the second 
kind. Our present paper can also be contrasted with recent work of Gelineau and Zeng [11] . who 
present an alternative approach to the combinatorics of the Jacobi-Stirling numbers. In a parallel 
development, Mongelli has recently established the total positivity of the Jacobi-Stirling numbers in 
[15] . In addition, in the recent manuscript [16] Mongelli shows that the Jacobi-Stirling numbers are 
specializations of the elementary and complete homogeneous symmetric functions; he also obtains 
combinatorial interpretations of a wide class of numbers which include the Jacobi-Stirling numbers 
as special cases. We remark that Comtet [6] further generalized the classical Stirling numbers in 
[6]. For an excellent account of Stirling numbers of the first and second kind, see Comtet's text [5j 
Chapter V]. 

The contents of this paper are as follows. In Section [2j we briefly review the Jacobi-Stirling 
numbers from the original context of left-definite theory. Section [3] deals with a comparison of 
various properties of the classical Stirling numbers of the second kind and the Jacobi-Stirling 
numbers. The proofs of many of these properties are similar to the proofs given in [3] so proofs, 
in this paper, will either be brief or omitted completely. In Section 0] we give a combinatorial 
interpretation of the Jacobi-Stirling numbers in terms of set partitions with a prescribed number 
of blocks. In Section [5] we study the Jacobi-Stirling numbers of the first kind and prove several 
properties of these numbers which are analogues of properties of the classical Stirling numbers of the 
first kind. In particular, we prove a reciprocity result connecting the Jacobi-Stirling numbers and 
the Jacobi-Stirling numbers of the first kind. In Section[6]we give two combinatorial interpretations 
of the Jacobi-Stirling numbers of the first kind in terms of ordered pairs of permutations with 
prescribed numbers of cycles. 

Notation. It is clear from (|l.lj) that {"} is a function of a + fi + 1, rather than of a and /3 

individually. With this in mind, we might wish to follow Gelineau and Zeng in setting z = a + /3 + l. 

It turns out, however, that it is more natural for our investigations to view these numbers as a 

z + 1 

function of 7 where 7 = — - — . With this in mind, we write {"} to denote the Jacobi-Stirling 
number {"} ^, where 27 — I = a + @ + I. With this notational change, we see from (jl.ip that 

n \ . _ 1)r+j (gr + 2 7 - l)r(r + 2 7 - 1) [r(r + 2 7 - 1)]" 

J'JV ^ r! ( j _ r ) ! r(j + r + 27) 

Observe that the numbers (that is, 7 = 1) are precisely the Legendre-Stirling numbers. The 
following table lists some Jacobi-Stirling numbers for small values of n and j. 



j/n 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


j = 


1 

















3 = 1 





1 


2 7 


4j 2 


87^ 


I67 4 


J = 2 








1 


67 + 2 


287^ + 2O7 + 4 


I2O7 3 + 1367^ + 567 + 8 


J = 3 











1 


I27 + 8 


1007^ + 1407 + 52 


j = 4 














1 


2O7 + 20 


J = 5 
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Table 1: Jacobi-Stirling Numbers 



A similar table, with the Jacobi-Stirling numbers in terms of a and j3, can be found in [10J; also, 
in [TT], the authors have a table of Jacobi-Stirling numbers given in terms of z. We show below, in 




JACOBI-STIRLING NUMBERS 



3 



Theorem I3,ll fiii). that the Jacobi-Stirling numbers satisfy a certain triangular recurrence relation 
that allows for a fast computation of these numbers. 

As with the classical Stirling numbers there occur, in a natural way, corresponding (signless) 
Jacobi-Stirling numbers [J] of the first kind (see Section [5] below; see also [7], [H], and p2|). In 

view of this, we occasionally will call {"} Jacobi-Stirling numbers of the second kind. 

2. Background 

We give a brief account in this section on the origin of the Jacobi-Stirling numbers. They were 
discovered in a certain spectral study of the Jacobi differential expression (j 1 . 3 j) : for specific details 
see [TO] , 

When 

A = A r = r(r + a + /3 + 1) + k (r € No), 
one classic solution of the Jacobi differential equation (II. 2p is 



(2.1) y{x)=pW\x)= ^ + ^F(-r,l + a + /3 + r;l + a;i^) (r G No), 

where Pr a '^\x) is the Jacobi polynomial of degree r as defined in the classic text of Szego |17t 
Chapter IV]. 

The most natural setting for an analytic study of (jl.2p is the Hilbert function space 

L 2 ((-l,l);w a , (x)) := L^(-l,l), 

defined by 



(2.2) L 2 at$ (-1, 1) := /:(-!,!) 



/ is Lebesgue measurable and J \f\ 2 w a ,pdx < oo 



with inner product 

-l 

-2 



(2-3) (f,9)a,(l ■= J i f(x)g(x)w a ^(x)dx (f,g € 1)). 

Pr > form a complete orthogonal sequence in L 2 o(— 1,1). 
Furthermore, in this space (called the right-definite spectral setting), there is a self-adjoint op- 
erator A^ a ^ in L^J— 1,1), generated by ^ Qj/ g[-], that has the Jacobi polynomials {-Pr "'^ }>?=o as 
eigenfunctions. Special cases of these polynomials include the Legendre polynomials (a = f3 = 0), 
the Tchebycheff polynomials of the first kind (a = (3 = —1/2), the Tchebycheff polynomials of the 
second kind (a = [3 = 1/2), and the ultraspherical or Gegenbauer polynomials (a = 0). 

The operator A^ a '^ is an unbounded operator but it is bounded below by kl, where I denotes 
the identity operator, in l? a J— 1, 1). Consequently, for k > 0, a general left-definite operator theory 
developed by Littlejohn and Wellman [14] applies. In particular, there is a continuum of Hilbert 
spaces {H[ a ' P) } t> Q where is called the t th left-definite Hilbert space associated with the 

pair {l? a o(— 1, 1), A^ '^). From the viewpoint of the general theory of orthogonal polynomials, it is 

remarkable that the Jacobi polynomials {P^" }^. f° rm a complete orthogonal set in H^ a '^ for 
each t > 0. The upshot of the left-definite analysis of (|1.3|) is that, for n € N, the integral composite 
power ^a[-] generates the n th left-definite inner product (-,-)^'^. In [ID], the authors prove the 
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following result which is the key prerequisite to establishing the left-definite theory of the Jacobi 
differential expression; it is in this result where the Jacobi-Stirling numbers are first introduced. 

Theorem 2.1. Let n G N. The n th composite power of the Jacobi differential expression (| 1 . 3j) . in 
Lagrangian symmetric form, is given by 



(2.4) r a>p [y](x) = _L-^(_l)i (4^(n,fc)(l -x)^(l + x)^+V j) (*)) ' (x E (-1,1)), 
where the coefficients Cj(n, k) (J = 0,1, ... n) are nonnegative and given by 

(25) c™(nk)-l° ifk = ° andc {l) (nk)-{ ^ * = ° 

In particular ) when k — 0, 

(2.6) = ( (!) (l - *r +J (l + s^W*)! (x e (-1, l)). 



7 



Furthermore, {"} is f/ie coefficient of x n 3 in i/ie Maclaurin series expansion of 



( 2 - 7 ) n i 



J l A l 



X < 



m(m + 27 — l)x V 7(7 + 27 — 1) 

m=l x 

From (|2.7p . we see that we can extend the definition of these numbers to include the initial 
conditions 

(2.8) jjjj =<5 n , and j°j = 5 jfi . 

The original motivation in the discovery of the Jacobi-Stirling numbers can be seen from the 
columns in Table 1. Indeed, the numbers in the n th column of Table 1 are precisely the coefficients 
of the n th power of the Jacobi differential expression £ a «[■]; for example, 

= -^W-l((l - x) a+3 (l + xf^)y'"{x))'" 
w a ,p\X) 

+ ((67 + 2)(1 - x) a+2 (l + xf +2 y"{x))" - (4 7 2 (1 - x) a+ \l + xf +1 )y'(x))'\. 

On the other hand, the rows in Table 1 demonstrate the horizontal generating function (12 .7p for 
the Jacobi-Stirling numbers. For example, reading along the row beginning with j = 2, we see that 

2 1 

TT 1 7 S 7TT = 1 + (67 + 2)* + (28 7 2 +20 7 + 4)t 2 + (1207 3 +1367 2 +56 7 + 8)i 3 + . . . . 

x — r(r + 27 — 1)1 

r=l x 

3. Jacobi Stirling Numbers Versus Stirling Numbers of the Second Kind 
The Jacobi-Stirling numbers {"} are similar in many ways to the classical Stirling numbers of 

the second kind, which we denote by {"}. As a first point of comparison, we note that, as reported 
in p3] (see also [8]), the Stirling numbers of the second kind are the coefficients of the integral 
powers of the second-order Laguerre differential expression m[-], defined by 

m[y](x) := {-x a+l e~ V(x))' (x E (0,oo)). 
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Indeed, for each n € N, 

(3.1) m n [y](x) = _L_^(_iy (l^j^e-^)^)) ' ! 

compare (|2.6[) and (|3.ip . 

The following table lists various properties of the Stirling numbers of the second kind. 



Property 


Stirling Numbers 2nd Kind 


Rational GF 


n l - rx ~ 2 \j\ X 

r=l n=0 w ' 


Vertical RR 


L"}-£C:i}'" 


Triangular RR 


{J} = {?:;}+ 

{0} = 5 n , ; {,} = <5j,o 


Horizontal GF 


x n = < . where 

3=0 ^ ' 
(x)j = x{x — 1) . . . (x — j + 1) 




Forward Differences 


A fc > for (n > j and fc € N ) 



Table 2: Properties of the classical Stirling numbers of the second kind 

Regarding the 'forward difference' entry in the above Table, recall that the forward difference of 
a sequence of numbers {x n }^ =0 is the sequence {Ax n }^ =0 defined by 

Ax n := x n+ i - x n (n G No). 

Higher-order forward differences are defined recursively by 

A k x n := A(A fc -V.) = ]T ( k \{-l) m x n+k - m . 

m=0 ^ ' 

Comtet [6j Propostion, p. 749] indicates the forward differences inequalities for the Stirling num- 
bers. Two excellent sources for the Stirling numbers of the second kind, and their properties, are 
the books by Abramowitz and Stegun jlj and Comtet [5, Chapter V]. 

We now state a theorem that has the comparable properties of the Jacobi-Stirling numbers; the 
reader will immediately observe the close similarities between Stirling numbers of the second kind 
and Jacobi-Stirling numbers. The details and proofs of most of these properties are given in [3] in 
the case of the Legendre-Stirling numbers (a = f3 = 0). Since the proofs are almost identical, we 
will only sketch proofs when necessary. 

Theorem 3.1. The Jacobi-Stirling numbers {"} have the following properties. 
(i) {Rational Generating Function) For all j 6 No, 

l = li_ r(r + 2 7 -l)x ^b/ 7 ^U7 7 V j(j + 2 1 -l))' 

in particular, for each n,j £ N, {"} is a polynomial in 7 with nonnegative integer coeffi- 
cients. 
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(ii) (Vertical Recurrence Relation) For all n,j £ No, 

(iii) ( Triangular Recurrence Relation) For all n, j £ N, 

n — 1 



{• 




f 71 — 








u- 














ft] 







+ i(i + 2 7 -i; 



and /or a// n, j € No, 

(n} cm 



ft}. 
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(iv) (Horizontal Generating Function) For all n £ No, 

n 

»• = £ 



j=Q 

v(7) 



(7) 



where ' is i/ie generalized falling factorial defined by 



w 3 



(7) ._ / 1 */i = 



Ilfci,(x-m(m + 27-l)) i/i € N. 
(v) (Forward Differences) For all k € No, 

A fe f U7 \ > ( > 

Proof. The rational generating function, given in (i), is discussed in Theorem 12.1) the complete 
proof of (i) is given in [10, Theorem 4.1]; we also refer to Section 4.2] where the authors 
determine this rational generating function in a different and simple way. The vertical recurrence 
relation in (ii) follows from (i); in fact, the proof is similar to that given in the Legendre-Stirling 
case which can be found in [3j Theorem 5.2]. Since the triangular recurrence relation is important 
in the combinatorial interpretation of the Jacobi-Stirling numbers, which we discuss in the next 
section, we give a proof of (iii). The initial conditions given in (iii) are part of the definition of 
{"} , given in fj2.8j> . From (i), we see that 



E i { J -! 1 }/- = a-^ +2 ,-iw|;{;[ 



n=3 

shifting the index on the sum on the left-hand side, and carrying out the multiplication on the 
right-side yields 

" " r n- j+1 



E{":i}^-i:{"} ^-a + 2 7 -i,e{"} 

n=j ^ J J 7 n=j ^ J J 7 n=j LJ J , 



n=3 1 n=j 



J 
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and this implies (in). The proof of (iv) is similar to the proof of the horizontal generating function 
for the Legendre-Stirling numbers given in O Theorem 5.4]. The proof of (v) is very similar to the 
proof given in [3j Theorem 5.1] in the case of the Legendre-Stirling numbers. □ 



For a different approach to parts (i), (iii), and (iv) see also [TT] , 

From Table 2 and Theorem 13. 1| observe that the rational generating functions for the Stirling 
numbers of the second kind (which we again note are connected to powers of the Laguerre differential 
expression in (|3.ip ) and the Jacobi-Stirling numbers (associated with the powers of the Jacobi 
differential expression in (12. 6p ) involve, respectively, the coefficients r and r(r + a + (3 + 1) in 
the denominators of these products. Remarkably, and somewhat mysteriously, these coefficients 
are, respectively, the eigenvalues that produce the Laguerre and Jacobi polynomial solutions of 
degree r in, respectively, the Laguerre and Jacobi differential equations. Computing the integral 
composite powers of both the Laguerre and Jacobi differential equations is entirely algebraic, and 
one would not initially expect these calculations to involve spectral theory. Furthermore, each self- 
adjoint operator in I? a J— 1, 1) generated by the Jacobi differential expression l a ^[-\ has a discrete 
(eigenvalues) spectrum only. It is natural to ask why the horizontal generating function specifically 
involves the eigenvalues {r(r + a + /3 + 1)}^L associated with the operator A^ a '^ (from Section[2]) 

that has eigenfunctions {P^ a ' }£1 an< ^ no * ^ ne se * °f eigenvalues of another self-adjoint operator. 
It seems that there is an interesting connection here that deserves further attention. 

The last point that we wish to make in this section is a very interesting connection between 
Legendre-Stirling numbers {"} and the classical Stirling numbers Sn of the second kind. In [3j 
equation (3.1)], the authors prove that 



It follows then that 



"} d 2n 1 ^, ,,j2j 

m=0 



j\ x {dxYidyY (2j)\ f- n V ; \mj 



p (j-m)x+{j+l-m)y | 

c \x=y=0 



1 f)2n 

- 1 u pJtt+GH-lWi _ p -(x+y)\2j I 

(2j)\ (dx) n (dy) n 1 ) lx=y=0 

d 2n ■ r- ^^ (e X+y - l) 2j 

(dx) n {dy) n (2j)\ lx=y=0 

Q2n 

= (&0W e "" e "°'~^ 

where (frj(-) is the vertical generating function for the Stirling numbers of the second kind; that is, 

~S« (expM-l)> 



n=0 J 
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Continuing, we find that 




Using the forward difference operator, this latter formula for the Legendre-Stirling numbers may 
be written in the compact form 

(3.2) {l} = = 

(3-3) h\r u ~ l] '' J '' xx: zrTp^ lx=y=0 when j - 2 ' 

with the obvious meaning of A x and Ay. It would be interesting to see if there is a similar connection 
between the Stirling numbers of the second kind and the Jacobi-Stirhng numbers. 

4. A Combinatorial Interpretation of the Jacobi-Stirling Numbers 

The Stirling number of the second kind {"} is the number of set partitions of {1, 2, . . . , n} into 
j nonempty blocks. That is, {"} is the number of ways of placing n objects into j non-empty, 
indistinguishable sets; for a full account of their properties, see [TJ and Comtet [SJ Chapter V]. 
With this in mind, it is natural to ask for a similar combinatorial interpretation of the Jacobi- 
Stirling number {"} . Indeed, Andrews and Littlejohn [2| have generalized the notion of a set 

partition to give a combinatorial interpretation of the Legendre-Stirling number {!•} ; one might 
call Andrews and Littlejohn's generalized set partitions Legendre-Stirling set partitions. More 
recently, Gelineau and Zeng [TTJ have found a statistic on Legendre-Stirling set partitions which 
allows them to interpret {") as a generating function over Legendre-Stirling set partitions. In 

this section we generalize Legendre-Stirling set partitions still further, to obtain objects we will call 
Jacobi-Stirling set partitions. When 7 = 1 it will be clear that the Jacobi-Stirling set partitions 
are in fact Legendre-Stirling set partitions, and we will show that for any positive integer 7, the 
Jacobi-Stirling numbers count Jacobi-Stirling set partitions. 

To describe a combinatorial interpretation of the Jacobi-Stirling number {"} , let [n\2 denote 

the set {li, I2, 2i, 22, . . . , n±, 71-2}, which contains two copies of each positive integer between 1 and 
n; we may say that these are the integers {1,2,..., n} with two different colors. By convention [0]2 
is the empty set. 

Definition 4.1. For alln,j,^ € No, a Jacobi-Stirling set partition o/[n]2 into 7 zero blocks and j 
nonzero blocks is an ordinary set partition of [n]2 into j+j blocks for which the following conditions 
hold: 
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(1) 7 of our blocks, called the zero blocks, are distinguishable, but all other blocks are indistin- 
guishable. 

(2) The zero blocks may be empty, but all other blocks are nonempty. 

(3) The union of the zero blocks may not contain both copies of any number. 

(4) Each nonzero block contains both copies of the smallest number it contains, but does not 
contain both copies of any other number. 

Example 4.1. As we see in the table below, there are 20 Jacobi-Stirling set partitions of [3] 2 into 
7 = 3 zero blocks and j = 2 nonzero blocks. 



Zero Boxes 


Nonzero Boxes 


Zero Boxes 


Nonzero Boxes 


0,0,0 


{l 1 ,l 2 ,3i},{2 1 ,2 2 ,3 2 } 


{3i},0,0 


{li,l 2 ,3 2 },{2 1 ,2 2 } 


0,0,0 


{l 1 ,l 2 ,3i},{2 1 ,2 2 ,3 2 } 


0,0, {3 2 } 


{li,l 2 ,3 1 },{2 1 ,2 2 } 


0,0, {3i} 


{1 1 ,1 2 },{2 1 ,2 2 ,3 2 } 


0,{3 2 },0 


{li,l 2 ,3 1 },{2 1 ,2 2 } 


0,{3i},0 


{1 1 ,1 2 },{2 1 ,2 2 ,3 2 } 


{3 2 },0,0 


{li,l 2 ,3i},{2i,2 2 } 


{3i},0,0 


{li,l 2 },{2!,2 2 ,3 2 } 


{2i},0,0 


{li,l 2 ,2 2 },{3i,3 2 } 


0,0, {3 2 } 


{1 1 ,1 2 },{2 1 ,2 2) 3 1 } 


0,{2i},0 


{li,l 2> 2 2 },{3i,3 2 } 


0,{3 2 },0 


{±1,1 2 }, {2 lj 2 2 ,3 1 } 


0,0, {2i} 


{li,l 2 ,2 2 },{3i,3 2 } 


{3 2 },0,0 


{li, 1 2 }, {21,22,3! } 


{2 2 },0,0 


{li,l 2 ,2i},{3i,3 2 } 


0,0, {3i} 


{1 1 ,1 2 ,3 2 },{2 1 ,2 2 } 


0,{2 2 },0 


{li,l 2 ,2i},{3i,3 2 } 


0,{3i},0 


{1 1 ,1 2 ,3 2 },{2 1 ,2 2 } 


0,0, {2 2 } 


{li,l 2) 2i},{3i,3 2 } 



As we show next, Jacobi-Stirling numbers count Jacobi-Stirling set partitions. 

Theorem 4.1. For alln,j, 7 G No, the Jacobi-Stirling number {"} is the number of Jacobi-Stirling 
set partitions of [n] 2 into 7 zero blocks and j nonzero blocks. 

Proof. For all n,j,j € No, let P(n, j, 7) denote the number of Jacobi-Stirling set partitions of [n] 2 

into 7 zero blocks and j nonzero blocks. Since the Jacobi-Stirling numbers are determined by the 

initial conditions and recurrence relation in Theorem 13. l( iii) . it is sufficient to show that P(n,j, 7) 

satisfies the same initial conditions and recurrence relation. 

Since the union of the zero blocks of a Jacobi-Stirling set partition cannot contain both li and 

1 2 , we see that P(n, 0, 7) = if n > 0. On the other hand, for each 7 there is one Jacobi-Stirling set 

partition of the empty set, so P(0,0,7) = 1. Finally, since the nonzero blocks of a Jacobi-Stirling 

set partition must be nonempty, we see that P(0, j, 7) = if j > 0. Therefore P(n,j, 7) satisfies 

the same initial conditions as | n ) . 

L 3 > 7 

To see that P(n, j, 7) satisfies the same recurrence relation, we note that for n > 1, Jacobi-Stirling 
set partitions of [n] 2 into 7 zero blocks and j nonzero blocks come in two disjoint types: 

(i) those in which n\ and n 2 are in the same box; 

(ii) those in which n± and n 2 are in different boxes. 

Each Jacobi-Stirling set partition in class (i) can be uniquely constructed from a Jacobi-Stirling 
set partition of [n — 1] 2 into 7 zero blocks and j — 1 nonzero blocks by appending a nonzero block 
containing only n\ and n 2 . Therefore there are P(n — 1, j — 1,7) partitions in class (i). On the 
other hand, each Jacobi-Stirling set partition in class (ii) can be uniquely constructed from a Jacobi- 
Stirling set partition of [n — 1] 2 into 7 zero blocks and j — 1 nonzero blocks by either inserting n\ 
into a zero box and inserting n 2 into a nonzero box, or by inserting n\ into an nonzero box and 
inserting n 2 into any box not containing n\. There are 7jP(n — 1 , j, 7) partitions of the first type, 
and j(7 + j — l)P{n — 7) partitions of the second type, so there are j(j + 27 — l)P(n — 1, j, 7) 
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partitions in class (ii). Therefore P(n,j,j) = P(n — l,j — 1,7) + j(J + 27 — l)P(n — 1, j',7), so 
P(n,j, 7) satisfies the same recurrence relation as {™} , and the result follows. □ 

Example 4.2. In [2l Example 4.4], the authors showed that {"}-, = 2 n ~ 1 . This argument general- 
izes to show that {"} = (27)™- 1 for all n G N. 

Example 4.3. In this example we give a direct combinatorial proof that { n ™ 1 }^ = 2Q) + 27Q), 
where (") denotes the usual binomial coefficient. By Theorem 14.11 the quantity {^x} is the 

number of Jacobi-Stirling partitions of [n]2 into 7 zero blocks and n — 1 nonzero blocks. In such a 
partition there must be exactly one number k which is not the minimum in its block, and at least 
one copy of k is in a nonzero block. If both copies of k are in nonzero blocks, then we may construct 
our partition uniquely by choosing the elements i < j < k of these blocks, and then placing k\ and 
&2 in blocks with minima i and j. Thus there are 2(g) of these partitions. Alternatively, if one 
copy of A; is in a zero block, then we may construct our partition uniquely by choosing the elements 
i < k of the nonzero block containing k, constructing that block with one copy of k, and placing 
the other copy of k in one of the 7 zero blocks. Therefore there are 27(2) of these partitions. 
Combining these two counts, we find that { n !^ 1 } = 2(3) + 27(2), as claimed. 

In a development independent of this work, Mongelli has recently given another combinatorial 
interpretation of the Jacobi-Stirling numbers [16J. In fact, Mongelli gives a combinatorial interpre- 
tation of a general family of numbers, which includes the Jacobi-Stirling numbers as a special case. 
Translated into our setting, Mongelli shows inductively that if z = 2j — 1 is a positive integer then 
{"} is the number of set partitions of \r1\2 into z zero blocks and j nonzero blocks for which the 
following conditions hold. 

(1) The zero blocks are distinguishable, but the nonzero blocks are indistinguishable. 

(2) The zero blocks may be empty, but the nonzero blocks are nonempty. 

(3) No zero block may contain the first copy of any number. 

(4) Each nonzero block contains both copies of the smallest number it contains. 

For convenience, we call one of Mongelli's set partitions a long Jacobi-Stirling set partition. 

For any n,j, 7 € No, we can give a bijection between the associated Jacobi-Stirling set partitions 
and the associated long Jacobi-Stirling set partitions. To do this, suppose a given long Jacobi- 
Stirling set partition has zero blocks Zq, Z\ , . . . , Z27-2 and nonzero blocks B±, . . . , Bj. To obtain a 
Jacobi-Stirling set partition, for each i, 1 < i < n, we do the following. 

• If i\ and ^2 are in different nonzero blocks, then we leave i\ and Z2 where they are. 

• If i\ and ^2 are in the same nonzero block and they are the smallest element of that block, 
then we leave i\ and Z2 where they are. 

• If i\ and ^2 are in the same nonzero block and they are not the smallest element of that 
block, then we put i\ in Z$ and we leave ii where it is. 

• If %2 € Z^ for < k < 7 — 1 then we leave i\ and 12 where they are. 

• If %2 £ Zry+k-i for 1 < k < 7 — 1 and i\ € B m then we put Z2 in B m and we put i\ in Z^. 

Now Z 7 , Z 1+ i, . . . , ^2 7 _2 are empty, so we remove them from our set partition. The result is 
a Jacobi-Stirling set partition, and we can reverse the process step-by-step to recover the original 
long Jacobi-Stirling set partition. 
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5. The Jacobi-Stirling Numbers of the First and Second Kinds 

We saw in Theorems 13.11 and 14.11 that the Jacobi-Stirling numbers are natural analogues of the 
Stirling numbers of the second kind (see also [11] and [IS]). For instance, both of these families 
of numbers have elegant horizontal generating functions: the Stirling numbers of the second kind 
satisfy 



(5.1) 



^ n = E 

j=0 



!:)• 



(n G N ), 



where (x)j is the falling factorial defined by 



(5.2) 



13 •- 



x(x - 1) • • • (x - j + 1) if j G N 
1 if J = 0, 



while the Jacobi-Stirling numbers satisfy 
(5.3) 



" = E 

3=0 



(") 



where (x)j is the generalized falling factorial defined by 



(5.4) 



x 



(7). = J1 if i=° 

3 ' \n£i(z-m(™ + 27-l)) ifjGN. 



In the case of the Stirling numbers we can invert (|5,ip to obtain the Stirling numbers of the first 
kind. In particular, these (signless) Stirling numbers of the first kind [™] may be defined by 



(5.5) 



£(-l)"+J 

j=0 



(n G N 



where (x) n is the falling factorial given in (|5.2p . Similarly, if we invert (|5.3p then we obtain a 
collection of numbers we will call the (signless) Jacobi-Stirling numbers of the first kind, and which 
we will denote by ["] . Specifically, 



(5.6) 



(x) 



(7) 
n 



j=0 



where (x)^ is the generalized falling factorial defined in (|5.4p . We immediately obtain the following 
biorthogonality relationships: 



m<j<n 
m<j<n 



n 




J. 





fn) 




"j" 




UJ 


t 


m 


7 



(n,m G N ), 
(n,m G N ). 



The table below lists the Jacobi-Stirling numbers of the first kind for small n and j; see also Table 
2 in Hi]. 
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j/n 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


j = 


1 

















i = i 





1 


27 


87 2 + 47 


487^ + 727 2 + 247 


3847 4 + 11527 il + IO567 2 + 2887 


3=2 








1 


67 + 2 


447 2 +527 + 12 


4007 tl + IOI67 2 + 7447 + 144 


3=3 











1 


I27 + 8 


1407 2 + 26O7 + 108 


3=4 














1 


2O7 + 20 


3=5 

















1 



Table 3: Jacobi-Stirling Numbers of the First Kind 

Like the Stirling numbers of the second kind, the signless Stirling numbers of the first kind satisfy 
a triangular recurrence relation (see [5j p. 214]): for all n,j G N we have 



n 




n — 1 


+ (n - 1) 


n — 1 






J ~ 1. 




J_ 




„ 3 . 



As we show next, the Jacobi-Stirling numbers of the first kind satisfy a similar triangular recurrence 
relation. 



Theorem 5.1. For all n,j € No we have 



(5.7) 

and for all n,j £ N we have 
(5.8) 







Vo 



and 



5 jV 



n 




n — 


1 


J. 


7 


J ~ 


1 



+ (n- l)(n + 2 7 -2) 



Proof. Line (|5.7p is immediate from (|5.6p . since (x)q 7 ^ = 1 and for all j G N the polynomial (a;}^ 
is divisible by x. 

To obtain (15.81) . first note that the coefficient of x J on the right side of (I5.6f) is (— lV +n T" - ] . Now 
observe that if n G N then 

(x)W = (x-(n-l)(n + 2 7 -2))(x>£ 1 

n-l 

= {x - (n - l)(n + 2 7 - 2)) ^(-l) n " 1+J 

3=0 

so the coefficient of x J on the left side of (15.61) is 



71—1 
J 



(7) 



-iy +n (n- l)(n + 2 7 -2) 



n — 1 
. j . 



+ (-l) J ' +n 



n — 1 

j 

n — 1 

i - 1. 



We obtain (I5.8P when we equate these two expressions for the coefficient of x J 



□ 



Next in this section, we prove a reciprocity result which connects Jacobi-Stirling numbers of the 
two kinds. To state this result, first observe that there is a unique collection [™] (n,j G Z) of 



polynomials in 7 satisfying the initial condition 



(5.9) 

and recurrence relation 
(5.10) 



^n,0; 



7 



n 




n — 


1 


j. 


7 


j ~ 


1 



+ (n-l)(n + 2 7 -2) 



n — 1 
J 



(n, k G 
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Moreover, these polynomials are the Jacobi-Stirling numbers of the first kind when n,j £ No- 
Similarly, there is a unique collection {"} (n,j £ Z) of polynomials in 7 satisfying the initial 

condition 

(5-11) {q} =<o, 

and recurrence relation 
(5.12) 



0j,0, 



+ i(i + 2 7 -r 



71—1 

j 



(n, k € 



Moreover, these polynomials are the Jacobi-Stirling numbers of the second kind when n, j £ No- 
It is not difficult to show that if n ^ and j / differ in sign then V]\ = {"") = 0, but we 

might hope that for various negative n and j we obtain interesting new polynomials [™] and {""} . 

Our reciprocity result, which is an analogue of a similar result |13t Line (2.4)] for classical Stirling 
numbers, shows that we actually obtain two familiar families of polynomials. 



Theorem 5.2. For all n,j £ Z we /iai>e 
(5.13) 



1-7 



Proof. The Jacobi-Stirling numbers of the second kind are uniquely determined by (|5.1ip and 
(I5.12p . so it is sufficient to show that the quantities L(n,j) = (— l) n +J [Z^\ also satisfy (15. lip 

and H532J). 

The fact that L(n,j) satisfies (15. lip is immediate from (15. 9p . so suppose n ^ and j ^ 0. Then 
we have 



L(n- 1J-1) 



-j + 1 
-n + 1 



1-7 



+ (-ir^(-j)(-J + l + 2(l-7)-2) 



1-7 



—J 
-n + 1 



1-7 



£(n,i)-i(j + 27-l)L(n-l,j), 



and the result follows. 



□ 



Concluding this section, we now turn our attention to the unimodality of the Jacobi-Stirling 
numbers of the first and second kinds. Recall that a sequence of real numbers {x n }^ =0 is unimodal 
if there exists two integers Aq < N2 such that (i) k < N\ — 2 implies < a^+i, (ii) Qa^-i < ajVi = 
• • • = cln 2 > OiV 2 +i) an( l (iii) ^ > -/V2 + 1 implies a& > a^+i- 

As in [H Section 5.7], we consider the horizontal generating function 



(5.14) 




(n € N ) 
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which, by Theorem 13,11 (hi) . is a polynomial of degree n with leading coefficient {"} = 1. Again, 
by this property, we find that Aq(x) = 1 and, for n > 1, 

^W=E{^Z 1 1 } ^' + £j0- + 27-l){ n ~ 
= x(A n ^(x) + 2 7 A' n _ 1 (x) + 

For example, we see that 

A 1 (x) = x, A 2 (x) = x(x + 27), A 3 (x) = x(x 2 + 2(37 + l)x + 47 2 ). 

Since 7 > 0, A\, A2, and ^3 are real polynomials, the zeros of which are all real, simple, and 
non-positive. Along the same lines as the proof of Theorem 5.7 in [3], we obtain 

Theorem 5.3. Let n E N. The zeros of A n , defined in Q5. 14\j , are real, simple and non-positive. 
Moreover, A n (ti) = 0. 

For a different proof, see [151 Theorem 4] . From Theorem 15. 3( equation (I5.6P , and a standard 
criterion for unimodality as given in Comtet [5} p. 270], we can state the following result. 

Theorem 5.4. The unsigned J acobi- Stirling numbers of the first kind and the J acobi- Stirling num- 
bers of the second kind are unimodal with either a peak or a plateau of 2 points. 

6. Two Combinatorial Interpretations of the Jacobi-Stirling Numbers of the 

First Kind 

The signless classical Stirling number of the first kind [™] is the number of permutations of 
{1,2, ... ,n} with j cycles. With this in mind, it is natural to ask for a similar combinatorial 
interpretation of the Jacobi-Stirling number of the first kind [ n ] . Indeed, Egge [7] has given a 

3 ■y 

combinatorial interpretation of the Legendre-Stirling number [™] in terms of pairs of permutations; 
one might call these Legendre-Stirling permutation pairs. More recently, Gelineau and Zeng |llj 
have found a statistic on Legendre-Stirling permutation pairs which allows them to interpret [ n ] as 

3 *y 

a generating function over these pairs. In this section we generalize Legendre-Stirling permutation 
pairs still further, to obtain objects we will call Jacobi-Stirling permutation pairs. These objects will 
come in two flavors, balanced and unbalanced. When 7 = 1 it will be clear that both the balanced 
and the unbalanced Jacobi-Stirling permutation pairs are in fact Legendre-Stirling permutation 
pairs. In addition, we will show that for any positive integer 7, the Jacobi-Stirling numbers count 
both the balanced and the unbalanced Jacobi-Stirling permutation pairs. In connection with these 
pairs, it will be useful to note that the cycle maxima of a given permutation are the numbers which 
are largest in their cycles. For example, if tt = (4, 6, 1)(9, 2, 3)(7, 8) is a permutation in S\o, written 
in cycle notation, then its cycle maxima are 5, 6, 8, 9, and 10. 

Definition 6.1. Suppose n,7 € No- A balanced Jacobi-Stirling permutation pair of length n is an 
ordered pair (7^1,^2) with t\\ € S n+J and 1T2 € S n ^-\ for which the following conditions hold: 

(1) 7Ti has one more cycle than 1x2- 

(2) The cycle maxima of tt\ which are less than n + 7 are exactly the cycle maxima of 1^2 ■ 

(3) For each k which is not a cycle maximum, at least one of vri(fc) and ^(fc) is less than or 
equal to n. 
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In [7] Egge defines a Legendre-Stirling permutation pair of length n to be an ordered pair (77, 7^) 
with 77 € SVi+i and TT2 £ <SVh such that 77 has one more cycle than 7T2 and the cycle maxima of 
77 which are less than n + 1 are exactly the cycle maxima of 7T2- If (77,7:2) is an ordered pair of 
permutations with 77 £ S n+ \ and TT2 £ which satisfies conditions (1) and (2) of Definition 16. 1| 
then the fact that 1:2 £ S n implies (vri,7r2) also satisfies condition (3). Therefore, the Legendre- 
Stirling permutation pairs of length n are exactly the balanced Jacobi-Stirling permutation pairs of 
length n with 7=1. The Legendre-Stirling permutation pairs are counted by the Legendre-Stirling 
numbers of the first kind ["] ; as we show next, the balanced Jacobi-Stirling permutation pairs are 
counted by the Jacobi-Stirling numbers of the first kind. 

Theorem 6.1. For all n,j,j £ No, the number of balanced Jacobi-Stirling permutation pairs 
(7Ti, 7T2) of length n in which 77 has exactly 7 + j cycles is ["] . 

Proof. Let a n j denote the number of balanced Jacobi-Stirling permutation pairs (77,7^) of length 
n in which 77 has exactly 7 + j cycles. 

It follows from (15.9P and fj5 . lOf) that [*] = 8^\. On the other hand, if (77,7^) is a balanced 

Jacobi-Stirling permutation pair of length 1 in which 717 has exactly 7 + j cycles, then we must 
have j < 1, since 77 £ #1+7 • If j = then some entry of 717 violates Definition 16. It 3) . so we must 
have j = 1. Moreover, when j = 1 both 77 and 7T2 must be the identity permutation, so aij = 5j \. 
Therefore the result holds for n = 1. 

Now suppose n > 1 and the result holds for n — 1; we argue by induction on n. To obtain 
a n ,ji first observe that by condition (2) of Definition 16.11 if (7Ti,7T2) is a balanced Jacobi-Stirling 
permutation pair of length n then 1 is a fixed point in 7Ti if and only if it is a fixed point in 7T2. 
Pairs (717, 7T2) in which 1 is a fixed point are in bijection with pairs (o~i, a 2) of length n — 1 in which 
o"i has j — 1 + 7 cycles by removing the 1 from each permutation and decreasing all other entries 
by 1. Each pair (77,7^) in which 1 is not a fixed point may be constructed uniquely by choosing a 
pair (01, a 2) of length n — 1 in which o~\ has j + 7 cycles, increasing each entry of each permutation 
by 1, and inserting 1 after an entry of each permutation. There are a n -ij pairs (01, 02), there are 
(n — l)(n + 7 — 2) ways to insert 1 so that 77(1) < n, and there are 7(71 — 1) ways to insert 1 so that 
77(1) > n. Combining these observations and using induction to eliminate a n -ij-i an d a n -\^ we 
find 



"n,3 



+ (n-l)(n + 2 7 -2) 



71 — 1 

3 



as desired. 



□ 



We have shown that the Jacobi-Stirling numbers of the first kind count balanced Jacobi-Stirling 
permutation pairs whenever 7 is an integer, but the orthogonal polynomials that give rise to these 
numbers include interesting special cases in which 7 is a half integer. Most notably, the Tchebycheff 
polynomials of the first kind occur when 7 = 5; while the Tchebycheff polynomials of the second 
kind occur when 7 = |. To address the combinatorics of the case in which 27 is an integer, we 
introduce unbalanced Jacobi-Stirling permutation pairs, and we show that these pairs are also 
counted by the Jacobi-Stirling numbers of the first kind. 

Definition 6.2. Suppose 71,27 £ No- ^ n unbalanced Jacobi-Stirling permutation pair of length n 
is an ordered pair (717,^2) with 77 £ S' n +27-i o,nd iT2 € S n for which the following hold. 
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(1) 77 has 27 — 1 more cycles than ix<i ■ 

(2) The cycle maxima of 77 which are less than n + 1 are exactly the cycle maxima of TT2 ■ 

It's not difficult to see that when 7 = 1 the unbalanced Jacobi-Stirling permutation pairs are 
exactly the Legendre-Stirling permutation pairs introduced by Egge, which are counted by the 
Legendre-Stirling numbers of the first kind ["] . As we show next, the unbalanced Jacobi-Stirling 
permutation pairs are counted by the Jacobi-Stirling numbers of the first kind. We note that this 
result is a special case of a result proved independently by Mongelli [I6j . 

Theorem 6.2. For all n,j,'y € No, the number of unbalanced Jacobi-Stirling permutation pairs 
(7Tx, tts) of length n in which 1^2 has exactly j cycles is [™] . 

Proof. Let a n j denote the number of unbalanced Jacobi-Stirling permutation pairs (717, 112) of length 
n in which 1x2 has exactly j cycles. 

It follows from and ([57L0]) that [V] = 6^7. On the other hand, if (717 , 7T2) is an unbalanced 

j 7 

Jacobi-Stirling permutation pair of length 1 in which 7T2 has exactly j cycles, then we must have 
j = 1, since 7T2 € S\. Moreover, when j = 1 both 77 and 7r2 must be the identity permutation, so 
a\ j = Sj %. Therefore the result holds for n = 1. 

Now suppose n > 1 and the result holds for n — 1; we argue by induction on n. To obtain 
a n j, first observe that by condition 2 of Definition 16.21 if (77,7^) is an unbalanced Jacobi-Stirling 
permutation pair of length n then 1 is a fixed point in 717 if and only if it is a fixed point in 7T2- 
Pairs (it 1, 7x2) in which 1 is a fixed point are in bijection with pairs (<7i, 02) of length n — 1 in which 
o"2 has j — 1 cycles by removing the 1 from each permutation and decreasing all other entries by 1. 
Each pair (77,7^) in which 1 is not a fixed point may be constructed uniquely by choosing a pair 
(o"l)0"2) °f length n — 1 in which 02 has j cycles, increasing each entry of each permutation by 1, 
and inserting 1 after an entry of each permutation. There are On-ij pairs (07,(72), an d there are 
(n — l)(n + 27 — 2) ways to insert our Is. Combining these observations and using induction to 
eliminate o re _ and a„_i^ we find 



+ (n-l)(n + 2 7 -2) 



n — 1 

j 



as desired. 



□ 



We conclude with an example involving unbalanced Jacobi-Stirling permutation pairs. 
Example 6.1. In this example we give a direct combinatorial proof that 



n— 1 



(n 



l)!ll( 2 7 + i) 

3=0 



By Theorem 16.21 the quantity ["] is the number of unbalanced Jacobi-Stirling permutation pairs 
(77,7^) of length n, where 77 £ 5 n+ 2 7 -i has 27 cycles and 7T2 G S n has 1 cycle. Since 7T2 has 1 
cycle, there are (n — 1)! choices for this permutation. Moreover, the cycle maxima for 77 must be 
n, n + 1, . . . , n + 27 — 1. Now we may place 1 in 7T2 in 27 ways, we may place 2 in 27 + 1 ways, 
and in general we may place j + 1 in 27 + j ways. Thus there are (n — 1)! 11^=0 (^7 + j) °f these 
unbalanced Jacobi-Stirling permutation pairs. 
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